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While the methodology of band structure unfolding has appeared in several publications, the 
original derivations of the unfolding formulas can be considerably simplified by using the fc-projection 
method. In this work, more transparent derivations of unfolded spectral weights are given by using 
the projection operator approach. A range of illustrative examples are also presented which include 
finite and random one-dimensional chains, and Kekule-textured honeycomb lattice. 
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I. INTRODUCTION 

The band structure property of energy spectra in a 
crystal depends on the possiblity of describing the sys¬ 
tem in terms of a primitive unit cell. In many elec¬ 
tronic structure calculations, however, it is necessary to 
use a supercell (SC) description, which results in folded, 
more densely packed and hard to interpret, bands in 
the smaller supercell Brillouin zone (SBZ). A challenging 
question is, therefore, how to recover, from the supercell 
calculation alone, the band structure in the larger primi¬ 
tive Brillouin zone (PBZ). An unfolding procedure must 
reproduce in the larger PBZ the original bands for a per¬ 
fect system, or give an effective band structure (EBS) for 
an imperfect one. The EBS has great interpretative value 
as it is comparable to the band structure of the perfect 
system and can reveal the perturbations caused by dis¬ 
order. The EBS is also important because it corresponds 
to angle-resolved photoemission spectroscopy (ARPES) 
measurements 

The methodology of band structure unfolding has been 
developed in a number of publications for different repre¬ 
sentations of the electronic structure.l^HSl The basic idea 
is to uncover the inherited symmetry of an SC eigen¬ 
function by comparing it with a reference Bloch func¬ 
tion of the underlying lattice. A key idea that facili¬ 
tates the conceptual description is to express the band 
structure in terms of the spectral function.l^El Suitable 
unfolded spectral-weight expressions that serve as un¬ 
folding formulas have been derived for plane-waves,the 
tight-binding method or Wannier functions ,1^ and linear 
combination of atomic orbitals (LCAO).I^ A few studies 
have circumvented tedious derivations of these formulas 
by defining t he s pectral weight formulas based on heuris¬ 
tic reasoning.!^ 

Various recent papers have applied group theory to 
band structure unfolding .ISHUl This approach, which has 
been referred to as the fc-projection method,li^ has the 
virtue that it does not require an artificial reference sys¬ 
tem, as it generates its own. Nevertheless, it leads to 
the same results as the previous approach based on the 
analysis of the spectral function.l^ 

Many other studies are using, and generalizing, the 
unfolding method to interpret the electron and phonon 


spectra of alloys, surfaces and disordered solids 

The purpose of this paper is to present simplified, more 
transparent, derivations of unfolding formulas based on 
the ^-projection operator. Following the theoretical de¬ 
scription in Sec. im a w ide selection of toy models are also 
presented in Sec. |III| that illustrate various applications 
of band structure unfolding. 

II. THEORY 
A. The projection operator 

A convenient notation is to use upper and lower case 
letters for SC and PC variables, respectively. A supercell 
is considered to contain a number AT of primitive unit cells 
at positions Yj. Correspondingly, a set of N reciprocal 
lattice vectors of the SC system, Gi, unfold a wavevector 
K in the SBZ onto wavevectors k = K -|- Gi in the PBZ. 
Using the translation operator, defined by T(r)^(x) = 
^(x -I- r), it was shown in Ref. [5] that the operator 

1 ^ 

Tk+g. = (1) 

J=i 

projects out of a supercell eigenfunction with wavevec¬ 
tor K a Bloch function with wavevector k = K -f G^. 
Equation 0 can also be obtained from the dehnition of 
projection operator in group theory, 

1 ^ 

( 2 ) 

where P is an irreducible representation of the group, 
which is for the translation group.^ The projec¬ 
tion operator is both Hermitian and idempotent. An 
elementary proof that Eq. 0 has the desired projective 
property is given in the following section. 

Spectral weight 

The norms of the projected components of an eigen¬ 
state, generally in the [0,1] interval, form the spectral 
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weights of the unfolding procedure, 

W"k,/(G,) = (TkV'KjIJ’kV'Kj) = (KJ|Tk|KJ). (3) 

B. Plane-wave representation 

To prove that Eq. Q is the required projection opera¬ 
tor, two conditions must be satisfied. First, the operator 
must project out of an SC wave function a PC Bloch 
function and, second, it must be idempotent. The proof 
below is based on a plane-wave representation. 

An SC eigenfunction can be represented as an expan¬ 
sion in plane waves, 

V’K(x)=^CK-cie*(*""°)•^ (4) 

G 

where the sum is over the set of all reciprocal lattice 
vectors of the SC system.l^ It is easy to verify the Fourier 
relation, 

1=1 g 

where {g} C {G} are the PC reciprocal lattice vectors. 
Applying the operator to wave function 0 , and using 
the Fourier relation ([^ gives 

1 Pk+g.i/’k(x) = ^ ( 6 ) 

g 

which is a PC Bloch function with wavevector K -|- G^. 
This proves the first condition. Applying T a second 
time, and observing that = 1, the second condi¬ 

tion, CP^ = CP, is also proved. This completes the proof of 
projective property. 

Equation Q can also be used in ([^ to ^e the spectral 
weights in the plane-wave representation,!^! 

t^KJ(G,) = ^|CK+G,-g^ (7) 

g 

The proof requires only the orthonormality relation, 

= ( 8 ) 

Within this representation, it is easy to verify the fol¬ 
lowing sum rule as well, 

ICK+G.-gl^ 

Gi Gi g 

= ^|Gk-g|^ = 1- (9) 

G 

Equation Q has a simple interpretation. The set {G} 
is partitioned into 3Nl sets, and the sums over each give 
the unfolded spectral weights. 


(a) (b) 





FIG. 1. Band structure of the ID chain in plane wave rep¬ 
resentation. The size of energy gaps is taken to be 5 units, 
(a) Reduced zone scheme, (b) Repeated zone scheme, (c) 
Extended zone scheme, (d) Unfolded bands. 

C. Atomic orbital representation 

In this section, unfolding formulas for localized ba¬ 
sis functions are derived. This covers both the lin¬ 
ear combination of atomic orbitals (LCAO), which uses 
non-orthogonal basis functions, and tight-binding models 
with orthogonal basis functions, i.e., Wannier functions. 
The LCAO method will be treated first as it covers the 
tight-binding method as a special case. 

The sites of the PC and SC systems are described by 
{r} and {R}, respectively. In the PC system atomic 
orbitals are denoted as |rm), with corresponding upper 
case letters used for the SC system. Furthermore, an 
orbital |RA/) in the SC is described as |R-f r(M), m{M)) 
in the PC. There are I and L unit cells in PC and SC 
descriptions, respectively, with I = 'NL. 

An eigenstate of the SC system is given in terms of an 
expansion in Bloch functions, 

|KJ)=^CK'^|K7V), (10) 

N 

where J is a band index and the Bloch functions are given 
by 

\^N) = ( 11 ) 

R 

The eigenvectors are solutions of a generalized eigen¬ 
value problem 

ijK ^ E^S{K)C^ (12) 

where the overlap matrix S'(K) is a Fourier sum of 
(0M|RA)=S'oM,RAr matrix elements reflecting the non¬ 
orthogonality of atomic orbitals. 
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Projection of an eigenstate (10) depends on the pro¬ 
jection of a Bloch function (11), which in turn depends 
on the projection of an SC orbital |RfV), 

1 ^ 

TkIRM) = -;^e-**^-'^-ir(r,)|R-fr(M),m(M)) 


3 = 1 
-N 


= -^e **"-'’^|R-l-r(M)-rj,m(M)). (13) 


3 = 1 

Since = e**"'*^, 

1 


Tk|KM) = —— +r(M),TO(M)) 

V R,rj 

= + r{M),Tn{M)) 


Having the projection of an eigenstate, Eq. (14), the gen¬ 
eral spectral weight formula, Eq. (Isl), can be applied to 
finds! 




KJ„ik-[r-r(M)] r. 
M ^ 


QN,rm{M)' ( 1 ^) 


NMr 


(a) (b) 




FIG. 2. Band structures of Peierls distorted one-dimensional 
chain, (a) Reduced-zone band structure from tight-binding 
model with the two hoppings t, t' = —0.9, —1.1. (b) Unfolding 
bands from same tight-binding model, (c) LCAO model with 
same hoppings as above, and overlaps s,s' = 0.15,0.25. (d) 
LCAO unfolded bands. 


A. Plane-waves in one dimension 


In the last step, the periodicity of SC was used to replace 
Sr ^ ^<5r.o- 

In orthogonal tight-binding, the overlap integral 
becomeJi^ 


^0N,rm{M) ^r(N)n{N),rm(M) 


= s, 


r(N),r ^n{N),m{M)') 


(16) 


so that the spectral weights, Eq. (15), simplify to 


^ NM 

X (17) 

For the LCAO and tight-binding models, in addition 
to Eq. ([^, there is another sum rule given bj!^ 

E^kj = ^, (18) 

.1 

where M is the number of orbitals in a supercell, so the 
right-hand side is just the average number of orbitals in 
a primitive cell. The proof can be understood by not¬ 
ing that a complete basis per PC must consist of M/AT 
orbitals. 


III. NUMERICAL EXAMPLES 

A number of toy models are presented in this section to 
illustrate the unfolding formulas of Sec.|n] The examples 
range from the simplest to illustrate the basic procedure 
to the more difficult ones of finite and random systems. 


As an illustration of Eq. 0. the simplest one¬ 
dimensional (ID) example is considered. 

The electrons move in a periodic ID potential of pe¬ 
riod a. The electronic structure is obtained by solv¬ 
ing the Schrodinger equation in momentum space (with 
h = 2m = 1 units), 


[{K - Gf - e.]CK-G + Y Ug'-gCk-g' = 0, (19) 


where Ug are the Fourier components of the potential, 
and G are integer multiples of 6 = 27r/a. For small po¬ 
tential the energy spectrum is free-electron-like with gaps 
of size 2\Ug\^ 

The band structure, shown in Fig. has three equal 
band gaps, which is obtained by setting Ub = U 2 b = 
Usb = U. The calculation is carried out by trun¬ 
cating the set of linear equations (19) to include only 


G = 0, ±&, ±26, ±35, ±4&, i.e., a 9 x 9 eigenvalue equa¬ 
tion. Figures [^a,b,c) show, respectively, reduced, re¬ 
peated and extended zone schemes as found in solid state 
textbooks.^ For this simple example the extended zone 
scheme is obtained by taking first band in first BZ, the 
second band in second BZ, and so on. However, this 
simple procedure is not feasible for more complicated 
examples in higher dimensions. Figure 1(d) shows the 
unfolded bands obtained by considering a hypothetical 
lattice with a period of a/4 and calculating the spectral 
weights according to (|^. Here the sum IC'fc-gP is ap¬ 
plied to the repeated zone calculation and includes only 
g = 0, ±45. In contrast to the abrupt changes in the ex¬ 
tended zone scheme, the unfolded bands are observed to 
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FIG. 3. Band structure of honeycomb lattice with Kekule 
distortion. The hoppings are t — —1.1 and t' — —0.8 produc¬ 
ing a band gap of 0.6 about zero energy. The special points 
r, K and M refer to center, corner and side of the hexag¬ 
onal PBZ, respectively. P and Q are special points of SBZ 
and correspond to K and M, respectively, (a) Reduced zone 
bands. Both K points of the primitive cell Brillouin zone are 
folded onto the T point of supercell Brillouin zone, (b) Un¬ 
folded bands showing the gap opening at the K point. Smaller 
gaps have also appeared at the zone boundaries of supercell 
Brillouin zone. 


diminish smoothly across the zone boundaries. 


FIG. 4. Energy spectrum of tight-binding electrons on a fi¬ 
nite chain, (a) Ghain of 11 atoms. Hoppings at both ends are 
enhanced by a factor of 2. Black points at K — 0 represent 
the eigenvalues, and the black curve is the exact dispersion of 
the infinite chain, E = 2tcos{ka). The K = 0 point has been 
unfolded to 50 equally spaced k points throughout the zone. 
The levels below and above the band limits each represent the 
overlap of two nearly degenerate eigenvalues, which can be re¬ 
solved on a finer energy scale (c). These are localized edge 
states and, thus, their unfolded points are extended through¬ 
out the zone, (b) Chain of 21 atoms. Increasing the number 
of atoms has made the unfolded spectrum to become more 
similar to the infinite chain band. Similary to the 11 atom 
chain, the edge levels are doubly degenerate on this scale. 


B. ID tight-binding model 


In this example, the ID chain is assumed to have un¬ 
dergone a Peierls distortion, which doubles the period to 
2a. The nearest-neighbor tight-binding model consists of 
two alternate hoppings t and t', and the transfer integral 
matrix is given by 


H^ = 


0 

t' + te 


2iKa 


t' + te 

0 


— 2iKa 


( 20 ) 


The two energy bands are shown in Fig. [^a). 

These bands can be unfolded to the BZ of the underly¬ 
ing cell of period a via the spectral weights 0- Writing 
the eigenvector coefficients C'^'^ simply as Cn the spec¬ 
tral weight of a band is given for this example by 

W = ^(C* + C'2*e*'=“)(C'i -b C'2e-*'^“). (21) 


integrals s and s'. This creates an asymmetry in the band 
structure. The overlap integral matrix to be used in the 
generalized eigenvalue equation is 


S(K) 


1 

s' + se 


2iKa 



( 22 ) 


The resulting bands are shown in Fig. [2 'c). 

The unfolded spectral weights, Eq. (15), can be written 
in a similar form as Eq. (21) as 


1 


W = ^(C^A + C;e"''°A*)(Ci + Cae 


—ika\ 


(23) 


where A = 1 + s e -p g' The unfolded bands are 
shown in Fig. [^d). 


The corresponding unfolded bands are shown in Fig.j^b). 


D. Honeycomb lattice with Kekule distortion 


C. ID LCAO model 

The LCAO model is similar to the tight-binding model 
but in addition to the hoppings t and t' there are overlap 


In the two-dimensional honeycomb lattice bond alter¬ 
nation is known as Kekule distortion, and has an impor¬ 
tant implication for the energy spectrum.!^ The struc¬ 
ture can be described by a commensurate lattice related 
to the underlying lattice by a nonsingular transformation 
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matrix with integer elements, 



The new unit cell is a (-s/S x '/3)R30° supercell with 3 
times the area of the primitive unit cell. (See Ref. [55] 


for more details.) There are two sites, or basis orbitals, 
per primitive unit cell and six sites per supercell. 

The Kekule distortion induces a gap in the spectrum 
of honeycomb lattice, given hy Eg = 2\t — t'\, where t and 
t' are the two hoppings. The transfer integral matrix is 
given by 




/ 0 t 

t 0 

0 t 

t 0 

0 

0 


0 

t 

0 

t'e"*K-A2 

0 

t 


The band structure, in reduced zone and unfolded ex¬ 
tended zone schemes, respectively, are shown in Fig. [^ 
This result is particularly graphic in showing the useful¬ 
ness of unfolding. While the SC folded bands have little 
resemblance to the bands of the perfect system, the un¬ 
folded effective band structure is quite comparable to the 
unperturbed one, different only in band gaps at the K 
point and the zone boundaries. 


E. Finite ID chain of atoms 

An important application of unf oldin g concerns the 
electronic structure of surface layers.l^l^ Commonly sur¬ 
face states of a crystal are studied by using a slab of 
several layers of atoms and forming a supercell by the 
addition of a vacuum layer, thick enough to prevent in¬ 
teraction between the slabs. For the vertical component, 
only the K = 0 point needs to be considered. 

The slab approach is illustrated here with a finite ID 
chain of atoms. The Hamiltonian matrix of a nearest- 
neighbor tight-binding model of electrons in a finite one¬ 
dimensional chain is given by 



t' 

0 . 

. 0 

0 


t' 

0 

t . 

. 0 

0 

0 

0 

t 

0 . 

. 0 

0 

0 

0 

0 

0 . 

. 0 

t 

0 

0 

0 

0 . 

. t 

0 


\0 

0 

0 . 

. 0 

t' 

0/ 


where the hoppings for the end atoms are allowed to be 
different. The eigenvalue for this Hamiltonian is solved 
to give the energy spectrum and eigenvectors, which cor¬ 
respond to AT = 0 in a supercell calculation. The results 
are shown for chains of 11 and 21 atoms in Fig. [^ The 
K — 0 point unfolds into a number Af (size of SC) of k 
points. The weights are seen to be smeared in fc as a re¬ 
sult of the chain being finite, but most of the weights are 


t 

0 

^/giK-A2 

0 

t 

0 


0 

Fe-*K-Ai 

0 

t 

0 

t 


t'e-*K-A3\ 

0 

t 

0 

t 

0 / 


(25) 


to 



FIG. 5. Unfolded band structures of a ID chain described by 
a supercell of 128 atoms. An ‘intrinsic’ broadening of 77 = 0.03 
has been used in order to obtain smooth spectral functions, 
(a) Perfect system with hopping t = — 1. (b) Disordered sys¬ 
tem with equal numbers of two hoppings given by t = —1.2 
and t' = —0.8. The spectral function was obtained by aver¬ 
aging over 12 random samples. 


concentrated near the dispersion curve of the infinite sys¬ 
tem. As the chain increases its length and approximates 
the infinite one better, spectral broadening in k is seen 
to decrease. There are two nearly degenerate states be¬ 
low, and two above, the band limits, which are smeared 
throughout the BZ. These are localized end point states. 


F. Disordered ID chain 


Another important application of unfolding is to obtain 
effective band structures of random alloys, which requi res 
very large supercells and is computationally expensive.^^^ 
This is illustrated here by the ID chain of atoms with two 
different hoppings distributed randomly. The approach 
used is to average the unfolded spectral functions of a 
set of random realizations of the system with a relatively 
small supercell. Alternatively, one can use one random 
realization of a very large supercell.^ The unfolded spec- 
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tral function for the ID model is given by 

A{k,E)=Y,WKjS{E-EKj), (27) 

,7 

where Wkj are the unfolded spectral weights given by 
Eq. ( [T7| ). In the numerical calculations, the delta func¬ 
tion is represented by a Lorentzian with width parameter 
77 . The calculated spectral functions of perfect and dis¬ 
ordered chain are shown in Fig. where the main effect 
of disorder is seen as a broadening in both k and E. 

IV. CONCLUSIONS 

Although band structure unfolding formulas can al¬ 
ready be used as a computational tool, it is desirable to 


have a clear and intuitive understanding of the concepts 
behind them. This revisit makes a contribution toward 
this aim by presenting a simple and unified description of 
the unfolding method based on the fc-projection operator. 
In addition, a range of simple examples are presented in 
pedagogical style that demonstrate interesting aspects of 
the unfolding method. 
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